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We have gone through a comparative study on two different kind of bulk viscosity expres¬ 
sions by using a common dynamical model. The Polyakov-Nambu-Jona-Lasinio (PNJL) 
model in the realm of mean-field approximation, including up to eight quark interactions 
for 2-1-1 flavor quark matter, is treated for this common dynamics. We have probed the 
numerical equivalence as well as discrepancy of two different expressions for bulk viscos¬ 
ity at vanishing quark chemical potential. Our estimation of bulk viscosity to entropy 
density ratio follows a decreasing trend with temperature, which is observed in most of 
the earlier investigations. We have also extended our estimation for finite values of quark 
chemical potential. 
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1. INTRODUCTION 

Ideal hydrodynamics successfully described the experimental data of transverse mo¬ 
mentum (pt) spectra and elliptic flow coefficient V 2 {pt) ® for different hadrons, 
produced in the heavy ion collisions. However, one can expect a weakly interacting 
gas like behavior of the medium, produced in the heavy ion collisions because of the 
asymptotic freedom of Quantum Chromo Dynamics (QCD). The argument went 
stronger when one coupled a 2 +i[2H11 

or dimensional ideal hydrodynamical 

model of Quark Gluon Plasma (QGP) phase to a hadron cascade one to prop¬ 
erly account for the viscous behavior in late hadronic stage.^i^i^ This eventually 
led to a concept of perfect fluidity^^ with an implication to be a strongly coupled 
plasma.fi^^^^ 
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In Au+Au and Cu+Cu collisions at Relativistic Heavy Ion Collider (RHIC), 
elliptic flow jg found to be large and not generally consistent with experi¬ 

mental observations upon consideration of different types of initial conditions. The 
absence of proper knowledge of initial conditions hints about the presence of vis¬ 
cosities in the QGP phase to reproduce V 2 data. It is observecP^^^^Mini that both 
shear and bulk viscosities suppress elliptic flow V 2 - However the smallness of this 
suppression leads to a prediction of very small specific shear viscosity valu^^*^ for 
the temperature region probed by RHIC or Large Hadron Collider (LHC). 

Bulk viscosity ( manifests itself by an addition of a diagonal term to stress 
tensor in the local rest frame. This happens due to local isotropic deviations 
from equilibrium. Bulk viscous pressure H becomes proportional to scalar expansion 
rate 9 at the location of fluid cell in the Navier-Stoke’s approximation, i.e. n(a;) = 
—(^d(x). So, expansion is opposed by the bulk term, manifested in the negative 
value. Thus, for an isotropically expanding fireball, bulk viscosity reduces the radial 
acceleration, thus decreasing radial flow. The picture is somewhat different in case 
of shear viscosity. Increase in radial flow by the shear term leads to flatter shape 
of particle px spectra,whereas the situation is reverse for bulk viscosity 
Alongside, in Ref.it has been found that for any reasonably small running coupling 
constant bulk viscosity is quantitatively less compared to shear viscosity. 

A list of references on the microscopic calculations of bulk viscosity coefficient 
is not very small.l^^lt^ Moore and SaremP^ and BucheP^ found that for systems 
with diverging specific heat near transition temperature Tc, the relaxation time 
can blow up with | still remaining finite. Coming to the expected nature of |, 
it should vanish for a system of massless non-interacting quanta due to conformal 
invariance. Presence of interactions in the system changes the scenario and leads to 
deviation from this conformal limit. However the deviation remains very small in 
all regions except near phase transition where strong interactions may lead to large 
correlation length.1^1^1^ Kinetic theory approaches estimate C to be varying as 
second power in this deviation. The technique has been applied in relaxation time 
approximatioiP^ as well as for systems of photons radiated by massive particles in 
equilibrium,!^ to which leading order QCD result^^also agree. For strongly coupled 
A/'=4SYM theory this variation however is taken to be linear in the deviation.!^ 
Computation for bulk viscous effects has been done by few authors considering 
hadron gaP^®^as well. Increase of | has been observed towards lower temperatures 
for massive pions,!^ in contrast to a decrease for massless pions.!^ However, near 
phase transition peak like behavior of specific bulk viscosity supports the general 
arguments,which happens because of long range correlations closely related 
to the chiral symmetry restoration. In presence of second order phase transition, | 
is supposed to diverge.!^!^ The critical behavior of bulk viscosity in Gross-Neveu 
and linear sigma models have also been explored in Ref.!^ and Ref.!^ respectively. 

In the Nambu-Jona-Lasinio (NJL) model, bulk viscosity coefficients C of quark 
matter has been calculated by Sasaki et al.,!^ Marty et al.,!^ Ghosh et al.!^ and 
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Deb et They have used the expression of bulk viscosity, which is based on 
the relaxation time approximation (RTA) in kinetic theory approacH^l^^^ or 
quasi-particle approach of Kubo formalism!^ Xiao et al.^ has calculated C in the 
Polyakov-Nambu-Jona-Lasinio (PNJL) model by using the expression based on the 
Kubo formalism with sum rule approach. Using the former and latter approaches, 
and Refs.l^l^^®^l^ have respectively estimated C in different 
other dynamical models. In this context, our interest lies in the study of both 
expressions of C under the framework of a common dynamical model and search 
some equivalence between two approaches. We have chosen Polyakov-Nambu-Jona- 
Lasinio (PNJL) model as our common dynamics to apply on both approaches of 
bulk viscosity calculations. 

The paper is organized as follows. Next section contains formalism part, where 
a brief discussion of PNJL model has been mentioned first and then the expression 
of bulk viscosity in two approaches are addressed. Next, in result section, the equiv¬ 
alence of two expressions of bulk viscosity is mainly investigated in numerical point 
of view for vanishing quark chemical potential (fi) and then results of finite fi have 
also been addressed. Finally we have made summary of our work. 


2. Formalism 


Kubo Formalisni^l^ relates viscosity coefficients to the correlation functions of the 
energy-momentum tensor. For the bulk viscosity coefficient ( of QCD matter, the 
relation is 


C = lim 

w ->0 


-ImG^(u;, 0) 
9a; 


lim 

o ;—>-0 


7rp(a;, 0) 

9a; 


( 1 ) 


where imaginary part of retarded Green function ImG^ and corresponding spectral 
density p, associated with C, are defined as 


-ImG^(a;, 0) = 7rp(a;, 0) = ^ dt j dVe*‘^‘([6»()(a:), 6»()(0)]) . (2) 

The structure of trace of energy-momentum stress tensor for QCD in low-energy 
theorems at finite temperature and density has been addressed in Ref.^^ Using 
standard Kramers-Kronig relation alongwith some basic thermodynamics and con¬ 
sidering the quark and gluon components, we geipSIHlIIl 


-G^ 


^(0,0) = 2 = -6( - flMl - flMl) + I6|e.| 3^ 


+ 





( 3 ) 


To extract bulk viscosity C, an ansatz is to be made for the spectral density p. 
Since the divergent contribution is deducted in the definition of quantities in right 
hand side of Eq.([3|) (realized in terms of the first term in right hand side), so the 
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high frequency perturbative continuurrP^ p{u) ~ is not included in left hand 
side of the same. In low frequency region, the following ansatz is chosen, 

0) ^ 9C ^0 

LO TT Wq + 

such that it satisfies Eq. O, where Wq is the mass scale corresponding to region of 
validity of perturbation theory. At wq >> T, it is expected that the spectral density 
will become perturbative and temperature-independent. So we obtain. 


c = 





+ 


rr. d d 


(mqq)^ + 6 -f -f 16|e„| 


( 5 ) 


which is the final expression for bulk viscosity coefficient, based on QCD sum 
rule. 

Next, we will come to another possible expression of which is based on the 
relaxation time approximation (RTA) in kinetic theory approach.^^^ Exactly same 
expression can also be derived from the quasi-particle approach of Kubo formalism, 
where spectral density of ^ can be expressed in terms of one-loop diagram.^^^®^!^ 
The idea of quasi-particle picture is introduced via inclusion of finite thermal width 
(r) in the internal lines of one-loop diagram. Instead of going to their detail deriva¬ 
tion, which one can see in Refs.l^^^ for RTA method and in Refs.^^^^ for quasi¬ 
particle approach of Kubo formalism, let us come to their final expression for For 
vanishing quark chemical potential (/r = 0), the RTA expression is 


c = 


12 

Y 


dpk h[l-h] 


( 2^)3 {Ek^T 


- c. 


-c: 




( 6 ) 


where the effect of Polyakov loop is absorbed into the PNJL distribution function 
designated by /$ .ESlEH ggi^ total bulk viscosity, we will have to take the sum 
over 2-I-I flavors. Here, = {k^ + is energy of quark and F = 1/r is its 

thermal width, which inversely determine its relaxation time r in the medium. 

Here in this work, we compute the bulk viscous coefficient using the method¬ 
ologies of LET and RTA. For our studies, we use the framework of the 2-I-I flavor 
PNJL model taking upto eight quark interaction terms. This model was developed 
by addition of Polyakov loop to the NJL model.^^^^^^ An insight into its formalism 
and recent developments can be found in Refs.^^^^^^^ 


3. Results and Discussions 

Let us start our discussion from numerical values of Q, obtained from the Eq. 
at /r = 0. These are plotted in the upper panel of Fig. O for two different values 
of Wo, which inversely control the numerical strength of C. Each of them exhibit 
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T (MeV) 



Fig. 1. (color online) as a function of temperature at vanishing chemical potential. Results of 
upper and lower panels are obtained from Eqs. m and m, which are based on sum rule and 
quasi-particle approaches of Kubo relation (RTA in Kinetic theory) respectively. 


visible peak structures near the transition temperature ^^169 MeV. The temper¬ 
ature dependence of quantities like e — 3P and (mqq)^ are collectively responsible 
for exploring this peak structure. The former quantity is closely related with speed 
of sound Cs- The deviation of from its limiting value, = 1/3, measures the 
violation of conformality, which is basically determined from the term (^ — cfj, 

associated with the former quantity. The second quantity (mqq'^rp vanishes in the 
massless limit and exposes a conformally symmetric medium. Hence, non-zero val¬ 
ues of both quantities make a link between the non-zero values of C and violation 
of conformality. The peak structure of C near indicates maximum breaking of 
conformal symmetry. 

Next, let us come to the results of C(T, = 0) from the Eq. ([6]), which are plotted 
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T (MeV) T (MeV) 

Fig. 2. (color online) Left: T dependence of masses m (lower panel) and (upper panel) of u or 
d (solid line) and s (dash line) quarks. Right: T dependence of normalized entropy density s/s^^ 
(upper panel) and square of speed of sound (lower panel), where SB limits of s (s = and 

(Cg = 1/3) are marked by red dotted lines. 


in the lower panel of Fig. O- Similar to Wq of Eq. ([S]), F = l/r in Eq. ([51) inversely 

controls the numerical strength of C- We have chosen two certain values of F or r, 

which can approximately reproduce the C, in same order of magnitude as done by 
Wo = 1 and 1.5 GeV in earlier case. On the basis of this approximate matching, we 
may get an equivalence between two parameters: 

T = 4 fm or, F = 0.049 GeV = wq = 1.5 GeV 

T = 7 fm or, F = 0.028 GeV = wq = 1 GeV , (7) 

which are coming from two different approaches. In this context, the reader should 
focus only on the peak structures near the transition temperature in both ap¬ 
proaches because their numerical strength (peak strength) are basically matched. 
Although, the numerical values of C, at temperatures away from Tc are quite different 
in two approaches. The possible reason is discussed in the next paragraph. 

One of the major difference between the two expressions of C is that the con¬ 
formal breaking terms are folded by (PNJL) distribution function of quark in the 
RTA expression (jH]) but it is absent in the LET expression ([S]). Now, this (PNJL) 
distribution function of quark always has a suppressing effect at low T after the 
folding or integration operation. So, this is the mathematical reason for lower val¬ 
ues of CiT) at T < Tc in RTA approach than that of sum rule approach. One may 
expect this kind of folding by the quark distribution function in the LET expression 
of C if one consider some different kind of ansatz for this spectral density, which 
may contain this kind of folding. In ref.it is distinctly mentioned that a possible 
uncertainty in this sum rule approach may be appeared from this ansatz for the 
spectral density. 

In the RTA approach, the terms of conformal symmetry breaking are (^ — c^) 
and {^fPrnQY Similar to earlier case, one can check again that both of the terms 
vanish in the massless limit. The role of '§r{'m-qq)j. in the sum rule approach is 
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1 1.2 1.4 1.6 1.8 2 


T/T^ 



Fig. 3. (color online) as a function of temperature at vanishing chemical potential in the 
sum rule (upper panel) and quasi-particle (lower panel) approaches. For comparison, Lattice-QCD 
values of ^ (squares with error bars) from Ref.^^ are pasted. 


equivalently played by the term, ^ = m? + Tm^, of the quasi-particle 

approach. To understand the temperature dependence of these quantities, we have 
plotted quark masses (lower panel) and their temperature derivative (upper panel) 
in the left-hand side of Fig. (ED- The peak structures of C are built by the peak struc¬ 
tures of which basically represent the transition between broken and restored 
phases of chiral symmetry. In the right hand side of Fig. the entropy density s 
(upper panel) and the square of speed of sound (lower panel) are plotted against 
temperature axis. Entropy density at high T domain reach to Stefan-Boltzmann 
limit(SB-limit), under 2-I-I flavor consideration. The s/s^g is plot¬ 

ted in the upper-right panel of Eig. m, where red dotted line is denoting the SB 
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1 1.2 1.4 1.6 1.8 2 


T/Tc 


Fig. 4. (color online) (upper panel) and ^ (lower panel) as a function of temperature for non¬ 
zero quark chemical potential. 


limit s/ssB = 1- Being inversely related with exhibits a dip near Tc, where 

rate of increment of s(T) is larger. Hence one naturally gets a peak structure in 
(i — c^) near T^, which contributes in the peak structure of SB limit of is 
marked by red dotted line in the figure. To display the approach of the thermody¬ 
namical quantities towards their SB limit, we have plotted Fig. [2] up to T = 400 
MeV to cover the high T region. 

Next, upper and lower panels of Fig. ([3|) show the temperature dependence of 
specific bulk viscosity C/s, obtained from two different approaches with the same 
parameters, taken for C{T) in Fig. ([T]). In both approaches, we observe decreasing 
nature of C/s{T) in the region of T > T^. Since C from sum rule approach is quite 
larger than that from quasi-particle approach at low temperature (T < Tc), so C/s 
increases to large numbers as one decreases the temperature at low temperature 
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Table 1. List of estimated values of f/s in earlier works at hadronic (2"^^ column) 
and quark (3'^'^ column) temperature domain. Dynamics with references and nature 
of temperature dependence are described in 1®* and 4**' columns respectively. 



VI 

Al 

Nature of T dependence 

LQCE)47I48mj 

- 

R 1 - 0 

Decreasing 

HRdSS 

^ 0.02 - 0.003 

- 

Decreasing 

HRG -1- 

Ri 0.02 - 0.1 

- 

Increasing 

HRdm 

Ri 0.01 - 0.035 

- 

Increasing 

HRG + 

Ri 0.025 - 0.12 

- 

Increasing 

hti]33 

- 

R 0.002 - 0.001 

Decreasing 

NJlIHl 

R^ 0.9-0.02 

R 0.02 - 0.002 

Decreasing 

NJimi 

R^ 1.7-0.13 

R 0.13 - 0.005 

Decreasing 

NJlJSSl 

- 

R 0.1 - 0.01 

Decreasing 

NJlHH 

R 0.61 - 0.11 

R 0.11 - 0.004 

Decreasing 

LSIVI^ 

R 0.61 - 0.11 

R 0.11 - 0.004 

Decreasing 

U nitarizatiorl^ 

R 0.04 - 0.027 

- 

Decreasing 

HRdMl 

R 0.15 - 0.025 

- 

Decreasing 


domain (T < T^). That is why the low T results have not been shown in the upper 
panel of Fig. (|31). In the lower panel of Fig. (|31), results of quasi-particle approach 
show peak structures at and hence, a discrepancy between the quantitative nature 
of C/s{T < Tc) in two approaches are observed. This is because the ({T < T^) from 
quasi-particle approach is quite smaller than that from sum rule approach. The C/s 
at T > Tc are more or less same in both approaches and quite satisfactorily agree 
with LQCD results, obtained in Ref.^ At high temperature domain, | in both 
approaches are tending to be zero. It indicates that both methodologies are obeying 
the general fact of high temperature QCD, which approaches towards the massless 
limit to attain its conformal symmetric behavior. 

We have extended our investigation for finite quark chemical potential by using 
Eq. ([5]) of sum rule approach in Kubo framework. In the upper panel of FigureQ, 
C is plotted for three different values of quark chemical potential fj,q= 50, 100, 150 
MeV keeping charge and strangeness chemical potentials {ffg and fis respectively) 
fixed at zero and Wq = IGeV. It is seen that with increase of fig, peaks of C 
are shifted towards lower T with higher numerical strength. Eor different values of 
chemical potential, the transitions take place at different temperatures in the QCD 
phase diagram, which causes shifts in the peak positions of bulk viscosity. 

The lower panel of Eig. (jH) shows the results for | as a function of T considering 
same set of values of chemical potentials, where we again see the increasing nature 
of I with increasing of fi. Similar trend was observed for shear viscosity to entropy 
density ratio in earlier studies,^ based on the same PNJL dynamics. 

To understand our results with respect to earlier works, we have made the ta¬ 
ble ll]), which contain approximate numerical values of //s, extracted from some 
earlier worksThe I*** column contains the information of 
dynamics for different works and next, the values of //s at hadronic and quark 
temperature domain for vanishing quark or baryon chemical potential are put in 
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2'^'^ and 3'’'^ columns respectively. The nature of the function, C/s(T) for different 
estimations are mentioned in the last column. We have also divided the table into 
two part, where references of upper part have used the expression of C, based on 
sum rule approach of Kubo formalism and the references of lower part have used the 
RTA or quasi-particle expression of Now most of the works observed decreasing 
nature of C/s(T) in both T < yjMlIlllIlHEol and T > yJlIlIlIlIlSIlssiET] j-egi^ie as well 
as in entire range of T jini[Ill[l41l56)[57] results, based on the Kubo expression of 
sum rule approach, are supporting this decreasing nature of C/s{T). In numerical 
point of view, one can relate our results to the divergence nature of ^/s near T^, as 
observed in Refs.^l^^®^ Now, on the other hand, the RTA expression is exhibiting 
a peak structure in the plot of Q/s vs T. This kind of peak structure is also ob¬ 
served in the Linear Sigma Model (LSM) calculation^SfllllZl as well as in NJL model 
calculations.^ We may assume indication of similar kind of peak structure from 
the increasing nature of (/s{T < Tc), observed in Refs.,^1^ whose calculations 
are based on the Hadron Resonance Gas (HRG) model withoulP^ and witll^l^ 
Hagedorn States (HS). So our investigations on two different expression of ( in 
PNJL model is revealing a qualitative discrepancy between the C/s{T < Tc) of two 
approaches, because of quantitative differences in their C(T < Tc). However, our 
C/s{T > Tc) are qualitatively very similar and for present set of values (ojq = 1, 1.5 
GeV and P = 0.049, 0.099 GeV) in both approaches match quite satisfactorily with 
results, obtained by Refs.^^^^ 

4. Summary 

In the framework of PNJL model with 2-1-1 flavor consideration, we have gone 
through a numerical investigation of two different expressions for bulk viscosity, 
based on two different approaches. One is Kubo formalism in sum rule approach 
and another is the quasi-particle approach of Kubo framework or the relaxation time 
approximation in kinetic theory approach. We have computed the temperature de¬ 
pendence of C and C/s for vanishing quark chemical potential in both approaches, 
where peak structure in C near Tc and the decreasing nature of C/s(r) (particularly 
at T > Tc) are their common numerical outcomes. The peak structure in C, which is 
interpreting the maximum breaking of conformal symmetry, is basically originated 
from two important quantities, associated with speed of sound and temperature 
dependent quark masses, which are adopted in both approaches. However, two ap¬ 
proaches contain two individual parameters, closely related with spectral width of 
Kubo formalism, which inversely control the peak strength of (. We have tried to 
observe their role for equivalence in C. At low temperature domain (T < Tc), we have 
noticed a numerical discrepancy in ((T) for these two approaches. The main reason 
is that the RTA expression contains the effect of folding by (PNJL) distribution 
function of quark, while it is absent in LET expression. 

Our finite density results of bulk viscosity are indirectly reflecting the fact of 
QCD phase diagram - the transition occurs at lower temperature as quark chemical 
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potential becomes non-zero and increases. This can be realized from the shifting of 
peak position in C, towards the lower temperature, when we increase quark chemical 
potential. 


Acknowledgments 

KS and SU would like to acknowledge Council of Scientific and Industrial Research 
(CSIR) and Department of Science and Technology (DST) for funding this work. SU 
would like to thank Avik Banerjee for useful suggestions. SG is financially supported 
from UGC Dr. D. S. Kothari Post Doctoral Fellowship under grant No. F.4-2/2006 
(BSR)/PH/I5-I6/0060. 


References 

1. J. Y. Ollitrault, Phys. Rev. D46, 229 (1992). 

2. P. Romatschke and U. Romatschke, Phys. Rev. Lett. 99, 172301 (2007). 

3. M. Luzum and P. Romatschke, Phys. Rev. C 78, 034915 (2008). 

4. H. Song and U. Heinz, Phys. Lett. B 658, 279 (2008). 

5. H. Song and U. Heinz, Phys. Rev. C 77, 064901 (2008). 

6. K. Dueling and D. Teaney, Phys. Rev. C 77, 034905 (2008). 

7. P. Bozek, Phys. Rev. C 81, 034909 (2010). 

8. A. K. Chaudhuri, J. Phys. G 37, 075011 (2010). 

9. B. Schenke, S. Jeon, and C. Gale, Phys. Rev. Lett. 106, 042301 (2011). 

10. P. Bozek, Phys. Rev. C 85, 034901 (2012). 

11. S. A. Bass et ah. Prog. Part. Nucl. Phys. 41, 255 (1998). 

12. M. Bleicher et ah, J. Phys. G 25, 1859 (1999). 

13. D. Teaney, J. Lauret, and E. V. Shuryak, Nucl. Phys. A 698,479 (2002). 

14. S. Ryu, S. Jeon, C. Gale, B. Schenke, and G. Young, Nucl. Phys. 904905, 389c (2013). 

15. H. Song, S. A. Bass, and U. Heinz, Phys. Rev. C 83, 024912 (2011). 

16. C. Shen and U. Heinz, Phys. Rev. C 83, 044909 (2011). 

17. U. Heinz and P. F. Kolb, Nucl. Phys. A 702, 269 (2002). 

18. M. Gyulassy, Structure and dynamics of elementary matter edited by W. Greiner 
et al., NATO science series H : Mathematics, Physics and Chemistry, Vol. 166 (Kluwer 
Academic, Dordrecht, 2004), p. 159-182 arXiv:nucl-th/0403032 . 

19. M Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 (2005). 

20. E. V. Shuryak, Nucl Phys. A 750, 64 (2005). 

21. H. Stocker, Nucl. Phys. A 750, 121 (2005). 

22. P. Bozek and I. Wyskiel, Phys. Rev. C 81, 054902 (2010). 

23. J. Steinheimer, J. Auvinen, H. Petersen, M. Bleicher, and H. Stocker, Phys. Rev. C 
89, 054913 (2014). 

24. U. Heinz, H. Song, and A. K. Chaudhuri, Phys. Rev. C 73, 034904 (2006). 

25. H. Song and U. Heinz, Phys. Rev. C 78, 024902 (2008). 

26. D. Molnar and P. Huovinen, J. Phys. G 35, 104125 (2008). 

27. E. Molnar, H. Niemi, and D. H. Rischke, Eur. Phys. J. C 65, 615 (2010). 

28. D. Teaney, Prog. Part. Nucl. Phys. 62, 451 (2009). 

29. U. Heinz, Relativistic Heavy Ion Physics, Landolt-Boernstein New Series, Vol. 1/23, 
edited by R. Stock (Springer Verlag, New York, 2010), Chap. 5. 

30. P. Romatschke, Int. J. Mod. Phys. E 19, 1 (2010). 

31. P. Kovtun, D. T. Son, and A. O. Starinets, Phys. Rev. Lett. 94, 111601 (2005). 








November 30, 2016 1:59 WSPC/INSTRUCTION FILE BulkTNJL 


12 K. Saha, S. Upadhaya, S. Ghosh 

32. G. Policastro, D. T. Son, and A. O. Starinets, Phys. Rev. Lett. 87, 081601 (2001). 

33. K. Dusling and T. Schafer, Phys. Rev. C 85, 044909 (2012). 

34. P. Romatschke, Eur. Phys. J. C 52, 203 (2007). 

35. H. Song and U. Heinz, Phys. Rev. C 81, 024905 (2010). 

36. G. D. Moore and O. Saremi, JHEP 0809, 015 (2008). 

37. A. Buchel, Phys. Lett. B 681, 200 (2009). 

38. A. Onuki, Phys. Rev. E 55, 403 (1997). 

39. K. Paech and S. Pratt, Phys. Rev. C 74, 014901 (2006). 

40. S. Gavin, Nucl. Phys. A 435,826 (1985). 

41. S. Weinberg, Astrophys. J 168, 175 (1971). 

42. P. Arnold, G. Dogan, and G. D. Moore, Phys. Rev. D 74, 085021 (2006). 

43. A. Buchel, Phys. Lett. B 663, 286 (2008). 

44. M. Prakash, M. Prakash, R. Venugopalan, and G. Welke, Phys. Rept. 227, 321 (1993). 

45. D. Davesne, Phys. Rev. C 53, 3069 (1996). 

46. J. W. Chen and J. Wang, Phys. Rev. C 79, 044913 (2009). 

47. D. Kharzeev and K. Tuchin, JHEP 0809, 093 (2008). 

48. F. Karsch, D. Kharzeev and K. Tuchin, Phys. Lett. B 663, 217 (2008). 

49. D. Fernandez-Fraile, Phys. Rev. D83, 065001 (2011). 

50. A. Dobado and J. M. Torres-Rincon, Phys. Rev. D 86, 074021 (2012). 

51. R. Marty, E. Bratkovskaya, W. Gassing, J. Aichelin, and H. Berrehrah, Phys. Rev. C 
88, 045204 (2013). 

52. Xiao Shi-Song, Guo Pan-Pan, Zhang Le, and Hon De-Fu, Chinese Phys. C 38, 054101 
(2014). 

53. C. Sasaki and K. Redlich, Phys. Rev. C 79, 055207 (2009). 

54. C. Sasaki and K. Redlich, Nucl.Phys. A 832 (2010) 62. 

55. S. Ghosh, T. C. Peixoto, V. Roy, F. E. Serna, G. Krein, Phys. Rev. C 93 (2016) 

045205. 

56. P. Deb, G. Kadam, H. Mishra, Phys. Rev. D 94 (2016) 094002. 

57. P. Chakraborty and J. 1. Kapusta, Phys. Rev. C 83, 014906 (2011). 

58. V. Chandra, Phys. Rev. D 84, 094025 (2011). 

59. V. Chandra, Phys. Rev. D 86, 114008 (2012). 

60. S. K. Das, J. Alam Phys.Rev. D 83 (2011) 114011. 

61. S. Mitra and S. Sarkar, Phys. Rev. D 87, 094026 (2013). 

62. S. Mitra, S. Gangopadhyaya, and S. Sarkar, Phys. Rev. D 91, 094012 (2015) 

63. G. P. Kadam, H. Mishra, Nucl. Phys. A 934 (2014) 133. 

64. G. P. Kadam, H. Mishra, Phys. Rev. C 92 (2015) 035203. 

65. G. P. Kadam, H. Mishra, Phys.Rev. C 93 (2016) 025205. 

66. G. Sarwar, S. Chatterjee, Jane Alam arXiv: 1512.06496[nucl-th]. 

67. J. Noronha-Hostler, J. Noronha and C. Greiner, Phys. Rev. Lett. 103, 172302 (2009). 

68. D. Fernandez-Fraile and A. Gomez Nicola, Phys. Rev. Lett. 102, 121601 (2009). 

69. D. Fernandez-Fraile and A. Gomez Nicola, Eur. Phys. J. C 62, 37 (2009). 

70. S. Ghosh, S. Chatterjee, B. Mohanty, Phys. Rev. C 94 (2016) 045208. 

71. H. B. Meyer, Phys. Rev. Lett. 100, 162001 (2008). 

72. R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 

73. R. Kubo, M. Yokota and S. Nakajima, J. Phys. Soc. Jpn. 12, 1203 (1957). 

74. 1. A. Shushpanov, J. 1. Kapusta, and P. J. Ellis, Phys. Rev. C 59, 5 (1999). 

75. H. Fujii and D. Kharzeev, Phys. Rev. D 60, 114039 (1999). 

76. P. Deb, A. Bhattacharyya, S. Dutta and S. K. Ghosh, Phys. Rev. C 79, 055208 (2009). 

77. S. K. Ghosh, S. Raha, R. Ray, K. Saha, and S. Upadhaya, Phys. Rev. D 91, 054005 
(2015). 








November 30, 2016 1:59 WSPC/INSTRUCTION FILE BulkTNJL 


A comparative study on two different expressions of Bulk viscosity in the PNJL model 13 

78. Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961). 

79. Y. Nambu, and G. Jona-Lasinio, Phys. Rev. 124, 246 (1961). 

80. T. Kunihiro and T. Hatsuda, Phys. Lett. B 206, 385 (1988). 

81. T. Hatsuda and T. Kunihiro, Phys. Rep. 247, 221 (1994). 

82. U. Vogl and W. Weise, Prog. Part. Nucl. Phys. 27, 195 (1991). 

83. S. P. Klevansky, Rev. Mod. Phys. 64, 649 (1992). 

84. M. Buballa, Phys. Rep. 407, 205 (2005). 

85. S. K. Ghosh, T. K. Mukherjee, M. G. Mustafa, and R. Ray, Phys. Rev. D 73, 114007 
(2006). 

86. S. Mukherjee, M. G. Mustafa, and R. Ray, Phys. Rev. D 75, 094015 (2007). 

87. C. Ratti, M. A. Thaler, and W. Weise, Phys. Rev. D 73, 014019 (2006). 

88. C. Ratti, S. Robner, and W. Weise, Phys. Lett. B 649, 57 (2007). 

89. W.-j. Fu, Y.-x. Liu, and Y.-L. Wu, Phys. Rev. D 81, 014028 (2010). 

90. W.-j Fu and Y.-l. Wu, Phys. Rev. D 82, 074013 (2010). 

91. A. Bhattacharyya, P. Deb, S. K. Ghosh, and R. Ray, Phys. Rev. D 82, 014021 (2010). 

92. A. Bhattacharyya, P. Deb, A. Lahiri, and R. Ray, Phys. Rev. D 82, 114028 (2010). 

93. A. Bhattacharyya, P. Deb, A. Lahiri, and R. Ray, Phys. Rev. D 83, 014011 (2011). 

94. A. Bhattacharyya, S. K. Ghosh, S. Majumder, and R. Ray, Phys. Rev. D 86, 096006 

( 2012 ). 

95. A. Bhattacharyya, S. K. Ghosh, A. Lahiri, S. Majumder, S. Raha, and R. Ray, Phys. 
Rev. C 89, 064905 (2014). 

96. A. Bhattacharyya, P. Deb, S. K. Ghosh, R. Ray, and S. Sur, Phys. Rev. D 87, 054009 
(2013). 

97. A. Bhattacharyya, R. Ray, and S. Sur, Phys. Rev. D 91 (R), 051501 (2015). 

98. A. Bhattacharyya, S. Das, S. K. Ghosh, S. Raha, R. Ray, K. Saha, S. Upadhaya, 

arXiv:1212.6010vl (2012). 

99. S. K. Ghosh, A. Lahiri, S. Majumder, M. G. Mustafa, S. Raha and R. Ray, Phys. Rev. 
D 90, 054030 (2014). 

100. S. K. Ghosh, T. K. Mukherjee, M. G. Mustafa and R. Ray, Phys. Rev. D 77, 094024 
(2008). 








